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Abstract: The quantitative structure-activity relationship analysis of a set of 31 peptide- mimetic analogues of HIV - 1 
protease inhibitors was performed by receptor surface analysis. RSA is a useful tool in situations when the 3D structure of 
the receptor is unknown since one can build a hypothetical model of the receptor site. A receptor surface model embodies 
essential information about the hypothetical receptor site as a three-dimensional surface with associated properties such as 
hydrophobicity, partial charge, electrostatic (ELE) potential, van der Waals (VDW) potential, and hydrogen bonding 
propensity. The surface points that organize as triangle meshes in the construction of the RSA store these properties as 
associated scalar values. Receptor surface models provide compact, quantitative descriptors, which capture three- 
dimensional information of interaction energies in terms of steric and electrostatic fields at each surface point, which in 
other techniques are calculated using probe interactions at various grid points. These descriptors can be used for 3D QSAR 
studies. 
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I. Introduction 

Molecules with different substituents are generated electronically with their biological activities. All the molecules 
were coded with a prefix "Exa". The receptor surface model is normally generated from the most active compounds in the 
data set. The rationale is that the most active molecules tend to explore the best spatial and electronic interactions with the 
receptor, while the least active do not and tend to have unfavorable steric and electrostatic interactions. Receptor models 
generated using the five best active compounds Exa4, Exa6, Exal2, Exal8 and Exa33 in the training set. After the receptor 
surface model has been generated, all the structures in the training and test sets can be evaluated against the model. The 
model can be used to calculate the energy associated with the binding of a molecule in the model. It can also be used to 
minimize a molecule by adjusting the geometry of the structure into a "best-fit configuration" based on the constraints 
imposed by the receptor model. This method creates hypothetical models, called receptor surface models; those characterize 
the active site of macromolecule based on the construction of surfaces to represent spatial and electrostatic properties of a 
receptor's active site. 

II. Result and Discussion 

In Receptor Surface Analysis (RSA), the major steps were (1) generating conformers and energy minimization; (2) 
aligning molecules using the MCSG method; (3) generating the receptor model; (4) evaluating the compounds in the 
generated receptor model; and (5) generation of equations by the genetic function approximation method. 

A receptor model was generated with the following: 

• Activity data 

• Soft receptor surface type 

• Energies calculated using electrostatic charge complementarity 

The option to minimizing molecules during evaluation was avoided. After the receptor model is generated, 
properties such as electrostatic potential and hydrophobicity were mapped on to the surface of the model. One property can 
be mapped at a time. Property maps were displayed as color regions on the receptor surface. 

These properties reflect the anticipated characteristics of the receptor that is being modeled. The intensity of color 
reflects the magnitude of the mapped property at a particular location. Properties that can be mapped include the following: 

a) Electrostatic potential: Each surface vertex is colored according to the potential value at the vertex position. Red areas 
have negative electrostatic potential, blue areas have positive potential and white areas have neutral potential. 

b) Charge: when this property is mapped, the surface color is based on the average of the charges of the template atoms 
closest to the receptor surface. Red areas are positively charged, blue are negatively charged and white areas neutral. 

c) Hydrogen Bonding: The color indicates the tendency for specific areas of the surface to act as hydrogen bond donors 
(green) or acceptors (blue). Areas of the model with no hydrogen bonding activity are colored cyan. 

d) Hydrophobicity: The surface is colored brown to map the hydrophobic areas of the model. Areas that are not 
hydrophobic relative to the scale on the panel are white. 
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Figure 1.1a shows the best active compounds embedded into the receptor surface model mapped with electrostatic 
potential; the red color represents negative energy values as favorable interaction sites, while the blue-colored regions 
represent positive energy values that are favorable sites for binding of the molecule on the receptor surface. Intermediate 
colors indicate the intensity and gray indicates neutral. Figure 1.1b shows the best active compounds embedded into the 
receptor surface model mapped with charge; the red areas are positively charged, the blue areas are negatively charged, and 
the white areas are neutral. Similarly, Figure 1.1c shows the best active compounds embedded into the receptor surface 
model mapped with hydrogen bonding; the blue areas act as hydrogen bond donors, the green areas act as hydrogen bond 
acceptors, and the pale blue do not have hydrogen bonding activity. Figure 1 .Id shows the receptor surface model mapped 
with hydrophobicity; the brown areas are hydrophobic. 

This color coding of the ligand-receptor interactions can offer a qualitative way of examining compounds, by 
introducing them into the virtual receptor and visually inspecting the favorable/unfavorable interactions; substituents that 
increase or decrease the binding affinity can be easily recognized, and one can make easily simple but accurate structure - 
activity estimations. Though the outcome is with statistically acceptable results for the training set, the prediction results of 
the test set are not equally encouraging. 




(Figure 1.1a: Electrostatic potential mapped on the receptor surface) 




(Figure 1.1b: Charge mapped on the receptor surface) 




(Figure 1.1c: Hydrogen bonding mapped on the receptor surface) 




(Figure 1 .Id: Hydrophobicity mapped on the receptor surface) 
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Four properties are independently mapped, generated and added them to study tables. Each study table was saved 
to do the QSAR model extraction. The relevant equations and their validation parameters are given here under. 

Electrostatic mapping: 

pIC50 = - 0.128124 - 1.33831 * "Inter VDW Energy" - 13.8893 * "VDW/693" - 27.1 182 * ELE/3708" 
Validation parameters are: 

R-squared = 0.884 

Cross validated R-squared = 0.799 

PRESS = 11.89 
Charge mapping: 

pIC50 = 4.19565 - 11.1954 * "VDW/645" - 32.9953 * "ELE/3838" - 16.9802 * "VDW/692" 
Validation parameters are: 

R-squared = 0.883 

Cross validated R-squared = 0.794 

PRESS = 12.22 
H-Bond mapping: 

pIC50 = 3.76507 - 9.89141 * VDW/644" - 19.1306* "VDW/692" - 33.215 * "ELE/2131" 
Validation parameters are: 

R-squared = 0.874 

Cross validated R-squared = 0.791 

PRESS = 12.36 
Hydrophobic mapping: 

pIC50 = 3.75074 - 17.9708 * "VDW/692" - 12.0473 * VDW/1545" - 38.0359 * "ELE/3706" 
Validation parameters are: 

R-squared = 0.885 

Cross validated R-squared = 0.814 

PRESS = 11.03 

The resultant models were of good statistical qualities having > 87 % explained variance and 79% predicted variance. 



Table 1.4: Predicted activity values based on RSA 



Code given 


Experimental 
pIC50 


Predicted pIC50 
from RSA 


Exal 


9.6 


9.17 


Exa3 


8.11 


9.25 


Exa4 


9.72 


9.7 


Exa5 


9.59 


9.67 


Exa6 


9.64 


9.54 


Exa7 


9.22 


9.64 


Exa8 


9.54 


9.23 


Exa9 


9.51 


9.48 


ExalO 


9.57 


9.25 


Exall 


5.53 


6.34 


Exal 2 


9.8 


9.27 


Exal 3 


7.56 


6.73 


Exal4 


9.14 


8.81 


Exal 5 


8.27 


8.89 


Exal 6 


9.28 


8.66 


Exal7 


9.6 


9.75 


Exal 8 


9.77 


9.94 


Exal9 


6.94 


6.75 


Exa20 


8.02 


8.4 


Exa21 


7.47 


8.07 


Exa22 


6.16 


6.77 


Exa23 


6.79 


6.57 


Exa24 


7.18 


8.93 


Exa25 


6.67 


9.19 


Exa30 


4.52 


4.24 


Exa31 


6.89 


6.45 


Exa32 


6.84 


6.69 


Exa33 


10 


10.22 


Exa34 


7.41 


6.71 


Exa49 


5.33 


5.4 


Exa50 


5.86 


6.34 
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III. Conclusion 

The descriptors VDW/693, ELE/3708, ELE/3838, VDW/692, VDW/644, VDW/692, ELE/2131, VDW/692, 
VDW/1545, ELE/3706, etc. are contributing at the grid point numbers in the space. The location and magnitude of a 
descriptor can be used as a guideline to improve the activity of molecules. 
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